Molecular ordering of precursor films during spreading of tiny liquid droplets 
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In this work we address a novel feature of spreading dynam- 
ics of tiny liquid droplets on solid surfaces, namely the case 
where the ends of the molecules feel different interactions to 
the surface. We consider a simple model of dimers and short 
chain-like molecules which cannot form chemical bonds with 
the surface. We study the spreading dynamics by Molecular 
Dynamics techniques. In particular, we examine the micro- 
scopic structure of the time-dependent precursor film and find 
that in some cases it can exhibit a high degree of local order. 
This order persists even for flexible chains. Our results sug- 
gest the possibility of extracting information about molecular 
interactions from the structure of the precursor film. 
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Since the discovery of dynamical layering by Heslot et 
al. there has been increasing interest in phenomena 
occuring at microscopic length-scales during spreading. 
The beautiful experiments of Refs. jl]-[I| reveal that both 
the molecular structure of the liquid and the type of sub- 
strate used influence density profiles of the droplets. For 
example, thickness profiles of tetrakis (2-ethylexoxy)- 
silane and polydimethylsiloxane (PDMS) droplets on a 
silicon wafer exhibit strikingly different shapes under 
spreading pi. Tetrakis exhibits clearly observable dy- 
namical layering, while the spreading of PDMS proceeds 
by a fast evoluting precursor layer of one molecular thick- 
ness. Furthermore, the experiments of Heslot et al. ||] 
and Valignat et al. Q also demonstrate the important 
role of liquid-surface interactions in determining the mor- 
phology of spreading droplets in microscopic scales. 

Because of these experiments theoretical work has 
gained new interest. Analytic theories to date deal with 
dynamical layering only A number of computer 

simulations have been performed [|s|-[l"5f to study the dy- 
namics of spreading with some emphasis on the molecular 
structure of the liquid. In particular, using Molecular Dy- 
namics (MD) simulations it was concluded in Refs. |l^ ] 
that the chain-like nature of the molecules can influence 
the structure of the precursor layer. 

The molecular structure of thin layers is of practical 
importance, too. One way of controlling the surface en- 
ergetics of a solid substrate is to use grafted molecules 
that adsorb on it, sometimes forming chemical bonds and 
brushed layers jlifl . Such surfaces have applications in 



e.g. coating and lubrication. Another important class of 
systems comprises amphiphilic molecules such as deter- 
gents where a strong asymmetry of interactions causes 
layered structures to form |17|. Spreading dynamics of 
such molecules is then of particular interest in trying to 
understand how these layers form, and how well ordered 
they are. 

Motivated by these considerations, we have chosen to 
study the dynamics of spreading of simple models of 
short, asymmetrically interacting molecules using MD 
simulations. This method is convenient since microscopic 
information about the structure and properties of the 
spreading layers is readily available in real time. Our 
model follows Ref. |12| where short chain-like molecules 
interact with each other via the Lennard-Jones (LJ) po- 
tential V(r) = 4e f ((a f / r) 12 — (af/r) 6 ), with parameters 
<7/ and ef for the width and depth of the potential, re- 
spectively. The n-mers comprise n LJ particles intercon- 
nected by a very rigid pair potential V c = (l/2)/c(r — tq) 2 , 
where k = 10000e//<7 2 .. 

The surface is modelled by a flat, continuum LJ ma- 
terial with a pairwise LJ interaction integrated over 
the volume z < to obtain an interaction potential 
Vi(z) = -A/z 3 + B/z 9 , where A = (2n/3)e s cr 6 s p, and 
B = (47r/45)e s cr s 12 p (T|, with p ~ 10 30 m~ 3 . Table 
I shows the parameters. Interaction asymmetry comes 
about through the choice of the potentials: Vi for one 
end of the chain (called the grafted end) and either Vu 
or Vui for the other. The potential Vi is much steeper 
and deeper than Vu and Vm, both of which are of equal 
depth. In addition, the equilibrium distance of the sur- 
face potential minimum as seen by the monomers is the 
same for Vj and Vu (the ordinary case), but for Vui (the 
shifted case) the minimum distance is about two bond 
lengths further away from the surface. For the shifted 
case this means that individual dimers tend to lie per- 
pendicular to the surface. 

The simulations are done at a constant temperature 
of kT s = 0.8e/, which is well above the triple point 
of an LJ material fll2]] . The dynamics is described by 
the equations of motion with a Nose-Hoover thermostat 
(|l2| , P~8 19 1 . They are solved using a modified velocity Ver- 
let algorithm ^0|. The Hamaker constant Ah is deter- 
mined by the effective bond length 6/ of the dimers; e.g. 
for bi w 100 A, A H « 10" 18 J which is a realistic value 
for surfaces [jl7| . Assuming that the effective mass of the 
molecules is about 10 5 amu Q, the time step in reduced 
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units (r.u.) is t r « 10~ 13 seconds. 

Fig. 1 shows snapshots of typical evolution of dimer 
droplets during spreading for the two cases. The random 
initial state of the ridge-shaped droplet is first carefully 
prepared, after which the surface interaction is turned 
on and spreading occurs perpendicular to the axis of the 
ridge. The dimers above the surface are mainly oriented 
perpendicular to it due to the stronger attraction of the 
grafted end. Dimers on the surface, however, behave 
quite differently. In the ordinary case the precursor film 
appears rather rough except very close to the edges of 
the droplet where the molecules have enough room to lie 
flat. In the shifted case the situation changes dramati- 
cally and the precursor film appears very well ordered at 
all stages of spreading. 

To quantify these observations, we have followed the 
time evolution of various density-density autocorrelation 
functions |20| . In particular, we have defined the precur- 
sor pair correlation function as follows: 

g(\f i -f j \,t)=N(t)- 2 <EE^)^)> • C 1 ) 

The total number of molecules within the precursor film 
is denoted by N(t), r~l(t) is the position vector of the 
centre-of-mass of the i th dimer, and the summations go 
over the dimers in the precursor layer. This function 
has been computed during the spreading dynamics. Af- 
ter a short inital transient, it assumes a characteristic 
steady-state shape for the times studied here. For the 
ordinary case, g has fairly sharp peaks corresponding to 
nearest neighbour dimers only which indicates that the 
layer is disordered and liquid-like. In the shifted case, 
however, clear and sharp peaks can be observed corre- 
sponding up to about fourth or fifth nearest neigbour 
dimers. The precursor layer in this case has a high de- 
gree of local order even at these elevated temperatures. 
We have also followed the evolution of an orientational 
distribution function which provides us with information 
about the relative orientations of neighbouring dimers. 
For the shifted case a sharp peak is observed correspond- 
ing to parallel orientations, which indicates that the ori- 
entations of neighbouring dimers are strongly correlated, 
too. In both cases the very late time regime corresponds 
to a diffusively thinning monolayer. 

Fig. 2 shows typical smoothed density profiles of 
dimers at three different times. In the ordinary case the 
profiles are fairly smooth and rounded, and no layering 
is present. At later times a step develops at the edges 
of the film where dimers tend to lie flat on the surface. 
In the shifted case the edge of the precursor film always 
remains very sharp and well-defined. It is conceivable 
that tendency towards layering on top of the precursor 
layer might be observed for larger droplets in this case. 

We have also studied the time-dependence of the pre- 
cursor film width w(t). Fig. 3 shows a comparison 
between the two cases. They are very similar, with 
a more rapid "almost linear" (~ i 9 ) region followed 



by a crossover to slower "diffusive" (~ t ' 5 ) behaviour. 
This is characteristic of the truly microscopic droplets 
studied here and has been observed in Refs. (H]|l^] as 
well. For our cylindrical geometry, within the first re- 
gion w(t) (x N(t) ex t - 9 , which for a true 3 — d droplet 
means that its radius R(t) ~ >/~N ~ t 0A5 since the flux 
into the first layer remains roughly constant. This overall 
R(t) ~ i 5 behaviour is in agreement with experiments 

m 

Furthermore, it is possible to estimate the correspond- 
ing effective diffusion coefficients in the two regimes. 
In the ordinary case for the "truly" diffusive late-time 
regime we find Di ~ 10~ 6 m 2 /s pl| . For the shifted 
case Di is about an order of magnitude smaller because 
the activation energy arising from neighbouring dimers 
is much larger due to the close-packed geometry of the 
layer. Most importantly, however, estimating a "diffu- 
sion" coefficient D e from the early-time (almost) linear 
region, we find that D e /Di w 100. This is in surprisingly 
good quantitative agreement with recent experiments, 
where the same ratio was measured between early-time 
and late-time "diffusion" coefficients |||. 

Additional studies of flexible chains consisting of four 
and eight particles where one of the ends is grafted with 
Vi reveal that the spreading becomes similar to the re- 
sults of Refs. Jl2| without any asymmetry of interaction, 
provided that the equilibrium orientation of an individual 
chain is parallel to the surface. In the case of a shifted po- 
tential, however, in which the equilibrium orientation of 
an individual chain is perpendicular to the surface, com- 
pact and well-ordered structure of the precursor layer 
still persists even for completely flexible chains. Com- 
plete results with different shifted potentials will be pub- 
lished elsewhere. 

Due to the microscopic nature of our droplets, it is dif- 
ficult to make comparisons with experiments. A spread- 
ing experiment with grafted polymer chains has been per- 
formed by Valignat et al. In some cases, the molecules 
become anchored on the surface and spreading stops. 
However, some of their density profiles do qualitatively 
resemble ours. For example, in the case of a microdroplet 
of PDMS-OH on a silicon wafer covered with a layer of 
behenic acid, a well-defined and sharp precursor film de- 
velops similar to our shifted case. We are not aware of 
any experiments probing microscopic order within the 
precursor film. However, it is interesting to note that in 
the experiments of Ref . f22|| , equilibrium ordering of thin 
layers of a simple liquid between two corrugated surfaces 
was observed. 

In conclusion, we have shown further evidence of how 
the surface-chain interactions can play a significant role 
in determining the microscopic structure of the precur- 
sor layer. In particular, through a simple model we have 
demonstrated how the internal structure of the film can 
depend on the asymmetry of surface interactions between 
the chain ends. This can lead to a high degree of local or- 
der if the molecules prefer a vertical equilibrium orienta- 
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tion. Additionally, our model reproduces the experimen- 
tally observed more rapid spreading at the early stages, 
which crosses over to slower diffusion at late times. Based 
on our results we suggest the possibility of obtaining 
information about the nature of these interactions by 
monitoring the density profile and precursor film during 
spreading. 
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TABLE I. Surface-monomer interaction parameters for 
the two cases discussed in the text. 



FIGURE CAPTIONS 



Fig. 1. Snapshots of a dimer droplet for the ordinary 
and the shifted cases. The grafted end is represented by 
a filled circle, (a) Initial configuration for both cases (as 
seen along the axis of the ridge), (b) The ordinary case at 
t = 30000 r.u., and (c) the shifted case at t = 30000 r.u. . 
Differences in the morphology of the precursor layer are 
clearly visible. 



Fig. 2. Smoothed density profiles for the two cases: (a) 
the ordinary case, and (b) the shifted case. The shoulders 
at the edge of the precursor film in (a) are due to dimcrs 
that fall flat on the surface. 



Fig. 3. Width of the precursor film w(t) vs. time 
for the two cases discussed in the text. This data is 
for 1525 dimers. Note the considerably slower late-time 
behaviour of the shifted case. 
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